GravitoMagnetic Field in Tensor- Vector-Scalar Theory 



Qasem Exirifard 

School of Physics, Institute for Research in Fundamental Sciences, Tehran 19538-33511, 7rar0 

We consider the perturbative regime of the TeVeS theory. We then compute the gravitomagnetic 
field that a slow moving mass distribution produces in its weak Newtonian regime. We report that 
the consistency of the TeVeS gravitomagnetic field near the Earth space-time geometry with the 
Lunar Laser Ranging measurements requires the scalar coupling -Jr- to be weaker than 3.3 x 10 -4 at 
the standard confidence level. This bound beside the astronomical data demands that the dimension 
full parameter of the theory, I, holds I > (1.65 ± 0.35) x 10 25 m. 
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The missing mass problem in the galaxies can be solved either by assuming the presence of dark matter, a form 
of matter yet not included and seen in the standard model of the elementary particles , or altering the newtonian 
dynamics and gravity. The ACDM model buys the first approach. The success of the ACDM theory has stemmed 
various searches to directly/indirectly observe dark matter particles and their signatures which can not be produced 
by altering newtonian dynamics/gravity. However no conclusive signals for the dark matter particles yet have been 
■ reported. 

\ The second approach started by the Milgrom theory of the modified newtonian dynamics (MOND) [l[ , then changed 
^vq to the aquadratic Lagrangian model for gravity 0] (AQUAL). This approach also has signatures that can not be 
produced by the dark matter paradigm, including the detection of the gravitomagnetic force in the MOND regime 
[J|. Now a generally covariant realization of the AQUAL model, the TeVeS theory Q, is known to represents the 
q-i modified gravity's approach. 

TeVeS equations of motion are more evolved compared to those of the Einstcin-Hilbert theory. Though its spherically 
^jPj 1 static solution Q - analog to the Schwarzschild solution- is known, the exact Kerr-like solution of the theory so far 
has not been reported. In one hand no solution of TeVeS possessing non-trivial gravitomagnetic field is known. In 
the other hand, we have acquired very precise empirical knowledge on the gravitomagnetic field around the space- 
time geometry of the earth. These include the Lunar Laser measurements Jy| and the gravity probe B Q. In this 
paper, we investigate how measuring the gravitomagnetic force near the earth constraints the TeVeS parameters. 
In so doing, we first review the field equations of the TeVeS theory. We next solve these equations at the linear 
£Nl . order around the space-time geometry of the earth. We show that the TeVeS gravitomagnetic force is less than the 
I/"") ■ Einstein-Hilbert gravitomagnetic field by the factor of |^ . Noting that the LLR measurements require the consistency 
with the Einstein gravitomagnetic field with the accuracy of 0.1% we bound the scalar coupling to be weaker than 
t— I ' 2.3 x 10 -4 at the standard confidence level. This value of k constraints the dimension-full parameter of the TeVeS to 
^ — 1 , I > (1.65 ± 0.35) x 10 25 m. These are the first constraints on I and k parameters of the theory. 

I. EQUATIONS OF MOTION OF TEVES 

X 

The dynamical fields in the Tensor Vector Scalar (TeVeS) theory [3] are the Einstein metric g ab , a time-like vector 
A a , and a scalar <j). There is also a non-dynamical scalar field a. The Einstein metric couples to the physical energy 
momentum tensor represented by T ab . The physical metric is defined by 

g ab = e~^g ah - 2A a A b sinh(2</>) (1) 

The vector field is forced to be the unite timelike vector with respect to the Einstein metric 

~g ab A a A b = 1 (2) 

We assume that matter is an ideal fluid. So the physical energy momentum tensor takes the form of 

Tab = pu a u b + p{g ab + u a u b ) (3) 

where p is the proper energy, p is the pressure and u a is the 4- velocity, all three expressed in the physical metric. 
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The metric equation of motion reads 

G ab = 8rrG(T ab + (1 - e-^)A^ l{a A b) + r ab ) + O ab (4) 

where G ab is the Einstein tensor constructed out from the metric and its derivatives, round bracket represent sym- 
metrization with respect to indices, e.g. A( a Bb) = A a B b + A b B a while r ab and G a fc stand for the contribution of the 
scalar and vector fields: 



Tab = <? 2 (<l>,a<f>,b - -g^h^gab - A^iA^b) - ^A"(f>, v g ab )j - hjl V F{kGa 2 )g ab , 



(5) 



e ab = K ( g^F^F* - -F^F^g ab ) - XA a A b (6) 



where F ab represent the field strength of the vector field 

F ab = d a A b - d b A a (7) 

and A is a Lagrange multiplier, K and k are respectively the vector and scalar couplings, I is a parameter and F 
stands for the TeVeS function. It is useful to define a function fx{y) by 

- /iFC/i) - \v 2 F'{v) = y (8) 
using which expresses the equations of the scalar fields to 

kGa 2 = pQsPhrhnhv) (9) 

{pOtPbr^J h af3 0, a ) = kG[g^ + (1 + e-^AM"]^ (10) 

where 

h ab = g ab _ A a A b ( n ) 

The equation of motion for the vector field reads 

KV b F ab + XA a + 8irG(T 2 A b V b (f> V> = 8ttG(1 - e-^)A b T ah (12) 

Recalling @, one can derive the equation of motion for the Lagrange multiplier by contracting the vector equation 
with A a . The A's equation is a linear algebraic equation in term of A. 

II. TEVES THEORY IN THE SPACE-TIME GEOMETRY AROUND THE EARTH 

The TeVeS theory coincides to the general relativity in two limits 

• Milgrom limit: K — > , I — > oo 

• Bekenstein limit: k — > 0, K oc k, I oc k~i 

The Bekenstein limit can be understood as the perturbative regime of the Milgrom limit. Here we consider the 
Bekenstein limit. We shall consider the Newtonian regime of TeVeS theory wherein fj, = 1 and the contribution of F 
to the equation of motion of the metric can be neglected Q • Note that the term neglected does not contribute to \ 
behavior of the fields around spherical mass distribution we are interested in. So we have 



T ab 



CT 2 (^>,a4>,b - ^g^^^^gab ~ (f> ^(A( a (j) fi) - -A" <p, v g a b)j (13) 

We are interested in the solution representing the space-time geometry around the earth. The linear equations of 
motion suffice to describe the dynamics of gravity near the earth with the precision required. So we ignore all the 
non-linear terms in the field equations. To this aim we consider that 



9ab — Vab + h ab 



(15) 
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and we keep the linear terms in h ab . We also utilize the fact that h ab oc T ab , <f> <C 1 and keep the linear leading terms: 

G ab = 8nGT ab -\A a A b (16) 
KV b F ab + XA a = (17) 
V^^-iMV^) = kG[g ab + 2A a A b ]T ab (18) 

Contracting ([TT]) with A a and recalling A a A a = — 1, results 

A = KA a V b F ab (19) 

inserting which back into (|17|) gives a second order non-linear equation for A: 

V^F aii + A»V V F^ A a = (20) 

We require spatial isotropy far away from the matter distribution. The isotropy demands 

A a = (1,0,0,0) (21) 
Af 0) = A a (22) 

The presence of the matter adds a correction to this value of A^ : 

A a (0) = A a + SA a (23) 

oo 

SA a = ^8 n A a (24) 

71=1 

where 8 n A a represents the term at the order of (h ab ) n . We are interested only in the linear correction to A : 8\A. 
Since the field strength of the A vanishes, 8\A a does not couple to h ab . Therefore 8\A a = which in turn implies 
that A is vanishing at the linear order approximation to the equations of motion: 

A = + O(h 2 ) (25) 

Due to the value of A (|25[) . the equation of motion for the metric (|16p simplifies to the ordinary Einstein equation: 

G ab = 8nGT ab (26) 

We are interested in the gravitomagnetism approximation to general relativity. So we consider the physical energy 
momentum of a pressure-less slow moving mass distribution where the non-zero components of the energy momentum 
tensor read 

Too = P (27) 
Toi = Ji = pv t (28) 

where p is physical density and Vi represents the physical velocity. Defining the Einstein gravitoelectric and gravito- 
magnetic potential as usual 

.goo = 1 + 2$ (29a) 
go* = Ai (29b) 



These quantities are solved to 

$ = $jv (30a) 
A t = Af H (30b) 

where is the Newtonian field potential and Af H is the gravitomagnetic potential of the Einstein-Hilbert theory. 

In the absence of the mass distribution, we require the physical metric to coincide to the Einstein metric. This sets 
the cosmological value of the scalar field to zero. This choice corresponds to the cosmological definition of time where 
in the physical time and the Einstein time coincides to each other. In the presence of mattter 4> also gets a series 
expansion in term of the deviation from the Minkowski metric 



J2 6 ^ ( 31 ) 



n=l 
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where 6 n (j> represents the term at the order of (h a b) n . We assume that m <C 1 so the equation of motion of 8i<j> can 
be approximated to 

V i V i Si(f> = kGp (32) 

This implies that 

M = (33) 

47T 

where <&jv as in the previous case represent the Newtonian gravitational potential. Due to the way that the scalar 
appears in the physical metric, beside the fact that A has only a temporal component, the contribution of J to 8i<f> - 
which we have neglected - does not contribute to the go; components of the physical metric. 

The physical gravitomagnetic and gravitoelectric fields are defined by the metric that a test particle probes, or the 
metric that appears in the worldline action of a particle 



dTg ^^ (34) 



where r is an affain parameter. It is the physical metric that leads to the physical gravitomagnetic and gravitoelectric 
force. Inserting (|33|) , (|29|) and (|30| into the definition of the physical metric ([lj , the non- vanishing components of the 
physical metric at the gravitomagnetism approximation follows 

goo - 1-2(1 + ^)$" (35) 

47T 

g i = A t (36) 
Therefore the gravitomagnetic and gravitoelectric field of the TeVeS theory read 

3>TeVeS = (1 + (37) 
4-7T 

jTeVeS = A EH (38) 

The value of the Newtonian constant is identified by fitting 1 /r 2 function on the gravitoelectric acceleration around 
stationary mass distribution. So the reported value of G is in fact G t, s = G(l + Let the gravitomagnetic/electric 
field be expressed in term of the observed value of the gravitational Newton constant: 

$TeVeS = *1 (39) 

jTeVeS = (J _ (40) 
47T 

This means that the gravitomagnetic field of the TeVeS solution representing the space-time geometry around the 
earth is weaker than the Einstein-Hilbert theory's prediction by the factor of 1^-. 

Ref. Q considers the lunar laser ranging measurements and reports that the Einstein-Hilbert gravitomagnetic force 
is consistent with the observation with 0.1% accuracy. This bound requires 

— < 3.3 x 1(T 4 (41) 

47T 

at the standard confidence level. This bound is two orders of magnitude weaker than what is commonly expected for 
k in order to keep the TeVeS theory in its MOND regime consistent with the cosmological/astronomical observation. 
However, a priori contradicting this expectation does not diminish the possibility to retain the theory consistent with 
the nature. The Gravity Probe B Q through the Pugh-Schiff effect [1,11] provides independent measurements of the 
geodetic and frame-dragging effects at an accuracy of 0.28% and 19%, respectively, neither of which improves the 
LLR measurements' constraint on k. 

Now that we have a direct constraint on k, we can constraint the parameter of I in TeVeS theory. Noting that 
K < 10~ 2 , the MOND limit of the TeVeS theory for the choice of the F function of Q, gives to the following critical 
acceleration 



v3fcc 2 

a = (42) 



The MOND terminology requires [T(| 



a = (1.0 ±0.2) x lCT 10 ^. (43) 



Our constraint on A; then demands that 

l> (1.65 ± 0.35) x 10 25 m (44) 
at the standard confidence level. This is the first constraint on the dimension full parameter of the TeVeS theory. 
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